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Abstract 



In previous work on quantum phase transition in granular superconductors, 
where mean-field theory was used, an assumption was made that the order pa- 
rameter as a function of the mean field is a convex up function. Though this is not 
{•"•^ . always the case in phase transitions, this assumption must be verified, what is done 

! in this article. 

\^ . Quantum phase transition in granular superconductors, modeled as Josephson - junc- 

G^ ! tion arrays (JJA) have been a subject of extensive investigations. For a comprehensive 

review see a book by Simanek [Q]. The picture of the phenomenon is as follows. At 
the bulk transition temperature Tq the magnitude of the order parameter of each su- 
perconducting grain becomes nonzero. When the energy of Josephson coupling between 



5-H . 

Qh. grains Ej is less than Tq, thermal fluctuations cause the phases 6 of the order parame- 

^ ! ter on different grains to be uncorrelated until the temperature is lowered to T = Ej. 

When the charging energy Ec, i.e. the electrostatic energy of a Cooper-pair located on 
I a fixed grain is comparable to Ej, the zero-point fluctuations of 6 destroy the long-range 



' superconductive order even at zero temperature. 



A lot of work have been done to study the phase diagram of the granalular super- 
conductor in {T / Ec, Ej / Ec) plane [Q]. Mostly mean field theory has been used, though 
recently a renormalization-group study was employed 0. In previous work, using mean- 
field theory, an important assumption was made: that the order parameter < cos^ > 
is a convex up function of effective mean field. The aim of this article is to verify this 
assumption. 

Quantum JJA are described by a Hamiltonian |l|]: 

H = Y.Pik{2efhihk -J2EikCos{di - 4). (1) 

ik ik 

Here matrix elements Pik are Coulomb interactions and En- are the Josephson energies 
between the i and k grains, = —id/d6k is the excess Cooper-pair number operator, 
conjugate to the phase 9^ of the k grain. While Pik are nonzero for all pairs. En. are 



1 



nonzero only for nearest neighbors. In the following we will consider a periodic array 
of identical grains with identical Josephson junctions between them. This idealization is 
reasonable when the disorder in Josephson energies and Coulomb interactions is not to 
big. In opposite case one faces a percolation problem. 

The mean-field Hamiltonian in Hartree approximation is obtained by replacing all 
operators, except two conjugate operators corresponding to chosen grain, by their average 
values: 

Hj^p = P{2efn^ - zfxEj cos 9. (2) 

Here z is the lattice coordination number; P = Pa is the diagonal element of the potential 
matrix; /i is the quantum-statistical and spatial average value of cos 6; the average value 
of fi is zero because of electroneutrality of a sample, and this is why the mean-field 
Hamiltonian depends on diagonal elements of potential matrix only. 
By introducing variables 

Ec = P{2ey;q = zfiEj/Ec, (3) 

Eq.2 may be rewritten as 

Hmf = -Ecid^de^ + qcose). (4) 

We denote the eigenf unctions of the Hamiltonian in Eq.(4) as \I',^ and \E'^ and its eigen- 
values as E^ and E^. We classify the eigenfunctions by condition "^ni^q = 0) = cos(ra6'), 
= 0) = sin(n6'). 

For theory to be self-consistent the value of the order parameter, obtained using the 
mean- field Hamiltonian should be equal to /i: 

_ En,aexp(-E°/T) <mcose\^'^> 
^ J:n,aexp{-E-/T) • ^ ^ 

From the Feynman-Hellman theorem one obtains 

<'^n\cos9\'^n>= -^dEjdq. (6) 

Ec 

This means that we just need to calculate E^^ and E^ to solve the problem. Hamiltonian 
in Eq. (4) is similar to the operator of the Mathieu equation : 

Huathieu = -d^/dx^ + 2q cos(2x) . (7) 
The relation between the eigen-values of (4) and (7) is: 

El^iq) = Ec\a2n{2qy, El{q) = Ec\hn{2q). (8) 

In Fig. 1 the order parameter /i, as a function of the mean-field q, calculated using Eqs. 
(5), (6), (8) is plotted. The eigenvalues of Mathieu equation a^iq), bm{q) are calculated 
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using a routine from the IMSL library, m's up to ten have been taken into account. For 
Fig. 1(a) only even order Mathieu functions have been used. This corresponds to taking 
into account only 27r-periodic i.e. those with integer n. For Fig. 1(b) both even and 
odd order Mathieu functions have been used, corresponding to 47r-periodic eigenfunctions 
of (4). Inclusion of 47r-periodic eigenfunctions may serve as an approximation to the 
continuous spectrum, i.e. to taking into account all \l/„ with real n. The latter situation 
corresponds to the case when the 2tt periodicity is broken by coupling of the Josephson- 
junction to an environment (for example by normal shunt resistance), which permits a 
continuous change of the charge on the junction. 

In the previous work only slope of /i(g) for small q was calculated using either expansion 
of Mathieu functions for g <^ 1 [0], or perturbation series for (4) 0]. Then superconducting 
transition temperature was calculated from: 

dlJ'iq,T)/dq\g=o = Ec/zEj, (9) 

assuming that fi{q) is convex up. Fig.l proves that this assumption is correct for both 
cases when only 2tt and when also An periodic eigenfunctions are included. 

One can consider two types of transition. First is at fixed ratio Ej/Eq and variable 
temperature. One can express the order parameter /i through the mean field q using Eq.3: 
/i = q/{zEj/ Ec), and plot this line in Fig.l. The intersection of this line with different 
curves fi{q,T) gives the value of the order parameter for each temperature. The second 
transition is at fixed temperature and variable ratio Ej/ Ec- Here one should plot several 
lines = q/ {zEj/Ec) for different Ej/Ec, and find their intersection with the fixed curve 
IJ,{q,T). As all curves fi{q,T) in Fig.l are convex up both mentioned above transitions 
are continuous. 

In Fig. 2 a phase diagram of the granular superconductor, calculated using Eq. (9) 
is presented. In addition the curve for continuous spectrum from Ref.[^ is given. A 
common feature of phase diagrams, obtained with 47r-periodic eigen-functions and with 
continuous spectrum is a reentrant transition. Indications for such a transition were 
observed experimentally, for example, in Ref.^. 

In conclusion, I have verified the assumption, made in earlier work, that the order 
parameter of granular superconductor is a convex up function of the mean-field, or in 
other words that the transition is of the second order. 

I am grateful to V.K Ignatovich and J.M. Kosterlitz, for useful conversations and to E. 
Simanek for correspondence. This work was supported by National Science Foundation 
Grant No. DMR-9222812. 
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Figure 1: Order parameter as a function of the mean-field q = zfj^Ej/Ec- In (a) only 27r 
periodic eigenfunctions of the mean-field Hamiltonian were used. In (b) also Air periodic 
ones were included. 
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Figure 2: The mean-field phase diagram of a granular superconductor. Sohd hne is 
obtained taking into account only 27r-pcriodic eigenfunctions of (4). Dashed line obtained 
by including also 47r-periodic ones. Dotted line is from Ref.[4], where a continuum of 
states was used. 
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